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ABSTRACT 

The problem of potential steady subsonic flow for lifting 
surfaces is considered. This problem requires the solution of 
an integral equation relating the values of the potential 
discontinuity on the lifting surface and its wake to the values 
of the normal derivative of the potential which are known from 
the boundary conditions. The lifting surface is divided into 
small (quadrilateral hyperboloidal) surface elements, 2T/ , 
which are described in terms of the Cartesian components of the 
four corner points. The values of the potential discontinuity 
and the normal derivative of the potential are assumed to be 
constant within each element and equal to their values* at 
the centroids of the elements. This yields a set of~linear 
algebraic equations. Numerical results are in good agreement with 
existing ones. 
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SECTION I 

FORMULATION OF THE PROBLEM 


1.1 Introduction 

References 1 and 2 present a general theory for compressible 
unsteady potential aerodynamic flow around lifting bodies 
having arbitrary shapes and motions. Reference 3 presents a 
general numerical formulation for complex configurations in 
steady subsonic flow. Results are presented in Ref. 4. However, 
such a formulation is not applicable to zero-thickness configur- 
ations (lifting surfaces) . The present work introduces a 
formulation suitable for use with lifting surfaces. 

The distribution of the perturbation aerodynamic potential 
^ , around a body of arbitrary shape is given by the following 
integral expression 


4I( E <{ 


where 



— - — 
r fm 


(r)l*x 


( 1 . 1 ) 


E = 0 inside the body 
E - 1 outside the body 
E = 1/2 on the body 

( 1 . 2 ) 

£s a surface surrounding the body and its wake, and n. 
represents the normal to the surface. 

If the distance between the upper and lower sides of 
the surface goes to zero (.zero-thickness body) , one obtains a 
lifting surface formulation 


r* JS 

t/ 


(1.3) 
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where 


7j - ^ ~ 


(1.4) 


The subscript u. stands for upper and l stands for lower. 
Equation (1.3) shows that the potential can be represented 
in terms of doublets on the body and on the wake. On the wake, 
the value of D is constant along a streamline and equal to D 
at the trailing edge. 


1.2 -Discretization 

By dividing the lifting surface into small elements (see 
also Ref. 3) and applying the mean value theorem for Eg. (1.3), 
one obtains 

^=21 h / ol ( 7 ) (1.5) 

where are suitable mean values within the element, and the 

summation is performed over the elements of the lifting surface 
and of the wake, which is approximated by straight vortex lines 
starting at the lifting surface trailing edge. The perturbation 
velocity , ir= at the point P^ , is given by 

n = = 

w (1.6) 


where 
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«W// V - 

Zp k " Vk 

* ( 1.7) 

is the velocity created by the element 21^ . The normal de- 
rivative at the point P ^ of the surface is given by 


(££) nt-VoV-ET Dj ni-nt 
J PzPk 


V 


( 1.8 ) 

The boundary condition to be satisfied at L points (L is the 
number of lifting surface elements) is 


V *7? =(l + v)> 


n - o 


which, when combined with Eg. ( 1.8 ) becomes 


( 1.9 ) 


SI Ak k ~ Bh 


for the L unknown Dl . In Eg. (1..10) 


(1.10 ) 


/)k 4 Vh 4 * ft k 

where i is the unit vector in the direction of the x-axis. 

The contribution of the wake elements adds only to the row of 
lifting surface elements in contact with the trailing edge. 
Once Eg. (1,10 ) is solved, the velocity VJp can be evaluated 
through Eg, (1.6) using the same coefficient^ 
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SECTION II 

HYPERBOLOIDAL QUADRILATERAL ELEMENT 
2a 1- Introduction 

Reference 3 introduces a new type of surface element, the 
hyperboloidal quadrilateral element, a short description of 
which will be given here. Then, the gradient of Eq. (.1.3) 

(the integral is obtained in analytical form in Ref. 3) will be 
computed and further, the result will be put in a simple vector 
form. 

2>2 Surface Geometry with Hyperboloidal Quadrilateral Element 
Let the geometry of the element be described by the 
vector 

p-~ p(j\ f 2- ) 

( 2 . 1 ) ' 

t t -r- 1 

where $ and £ are the generalized curvilinear coordinates 
(Fig. 1) . 

The two base vectors are given by 

. -> 

-* OP 
CLC “ T 

^3 ( 2 . 2 ) 

and the unit normal to the surface is obtained as 

r* 

I 5i * ^a.1 (2.3) 

The surface element tft^Tis 

dz. - la’t 4? 1 X a.1 df j C |7, xa t \d f 'd^ 


(2.4) 
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The expression for r is 


l 


r x ' *k ' 
z - ?h 


(2.5) 


Now, consider the equation 

-> /~ 1 <I -- 1 

Pc ♦ ??;(.,*, tl , B .„ 

^ 

The above equation represents a hyperboloid (Fig, 2), Pc 
represents the centroid of the element 21 ^ , with "§ - f - 0. 

The corner points of the element , are P*^, P+-— , 

, P , and they are fed in as geometry inputs in the computer 

program implementing the theoretical formulation. The relation- 
ship between the corner points and P c , , P^ #■ P 5 f is 




r “1 


f ' 

Pc 


11 1 1 


P ++ 

p, 


11 -1 -! 


p + _ 

<-> 

1 - 


( 

-H - 

! P, 


1-1 1 -1 


P- + 

£ . 


1 -1 -1 i 


P — 


( 2 . 7 ) 


2.3' The' Doublet Integral 

Looking again at Eq. (1.13), it can be written in the. 

4 

following form 

- fv I) -fir ( - Iv'f-t, 1) 


where 


Zr flf)=-zn Vo I3>(L 1) 


( 2 . 8 ) 


(2.9) 
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where the doublet integral 1^ ( ^ ^ ) was obtained in analytical 

form {Ref. 3, Eq. 6.6) &S 
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talVfr-a, xaif _ t 

2 ■ £ (i ci i xa l )V i • 2 * ft i. 1 ) I llVi ■ &' K a 0 

( V x 5) • £* «j. + 2 x a, -1 

(txZt ■ %* %■ a,x ai)-+ V - a, xQi %■ f j- 


iSuatKxa^rui 


(l>* a, • £x4t ■+. £<«/ £•£ 1- a, sail - 


■ (<i < a ,■ ^x.aJ)( £-V <£- <ii x + V • x <*i £• £) 

( 2 . 12 } 

Next, it is convenient to introduce some classical concepts 
of tensor analysis. Consider the relationships between the two 
sets of conjugate base vectors 



— / & i x 0 . 3 

a - — ~z — ~ 
Cl t • * x # 3 

and 



-■_ . - <3 3 X & 1 . 

1 ci'-ct^xci 1 i z a* . a^xa 1 1 

f 

where a-^ « n. Note that 

-a 2 - a z -n ( 2 . 15 ) 

and 

1 

oi.XO-l j 




(2.16) 
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Furthermore , using classical notations/ it is possible to 


write 


9-V'a, + /a, + = P'a/ 


= i?, a' + ^ p. a 


(2.17) 


with 


V? = P • a' ' 


\) t - - V • CL( (2.18) 

Moreover/ it is convenient to consider the three derivatives 
Q an d Ql$/Or) . Since $ is an arbitrary 

vector, Eq. (2.12) with V =■ & > , yields 


O Id 


|$F K Gt\ a zl*~ Id 


(aixa/fxQi + ^xara^ai)^^ 2 - a, - 


a^. )(2-ai 2. -0.iX.ai_ + a, a, xal 1- 2 ) 




(t-Za'ra*- %■&, £• O-x) t a { 


T=4r'-^:i- t ,2. .— .- a,-a l -(<i-a / ) z ] 

a.xa,| 5 -|'t*4 i | It A L 


a • xa ?. 

141 


$ 

:r T' 1 2. 


(2.19) 
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Similarly, for - XX ^ 






- - g/.a.xa*- 


( 9L >*«•,• ^xQi)f 9.-a<xat-f ai-^^xa^ 2-^ 


SL • °- / x 


jcLxa,| i -|2xaj.i 2 -|e-J 


( «2 • £ 5 ,- a i - a L ■ a z % ■ a t ) q_ ■ % - 


( q_ <i <x r a z - <z-a,<i - a t ) a z - 


£ CLf K XL X 2. • a 


15.1 

Finally ,'\ for V -ft and using Eqs. 
yields, 


(2.20) 

(B.l) and (3.45) of Ref. 3 5 


O 1 7) 


— I > sr — t 2. \ •“ 


^ V 


ir K* a 'llt* a *! M\ = 


(n *.a, ■ £x<?, • n xa>_ ) a . j? %-a,xax 


( 9 . *■ a > • 1 * a *)( t i~ '< - «/ x «i -+■ n • a, a 2 ■ £ J — 

(%•?> a, -cit-n • «1 c?-<z, -+ j-n cC r a t ~f-a: z a r 7i)x 

- i’2 2- M«! - (T-f <?•«, Z-a-ijf-n %-a,xa 2 ~ 

■ ( i- 3 cTi • <u - 2 • <2, <j • a z ) n ■ a, x 5 "j f- f" = 

( 5- 2 «/ •«! + J-Q'2-ai)t-n %'a,xa L - 



‘ CK 
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i'&i 1 n-a/XGi — 


[( £• 2 6, •«! + 2 - a, Q'?j)( < 2' a, k a 


a , x a t. 


(“£* ICLt-a--*.- ?•«/ ?-«0i a '' a ^r^' X 1 


4 


s, *&u 


ja ,Aa i| - j 5 L a, *ai.a,x2-2-a,c3ix«<-A 1 * 4 ] 


+ fl-i 5 ,- a- L [ , l'?.|ai xaxl-2-ai a, taiami -%• a.i a^xa.,- ai*z]~ 


(V y_ a, -a z - <z- a, %. a*) %■{ }a, ■ &i} Z '] - 


) a, x a 


~ j 2. a , f-« v 2- 2 [*/ :«./ - /a, ■ 2F*) 1 ] - 


a , 9- <2 ^ . 4 a v fa r ^i ?■ ,. <3 ~t-g -(3 j (a±-Ai o- a ; -g ■ <s\ &,• a 2}] - 

\ 

--—-rr-j- J f 

^ [?‘ a < ai-ai-f-aia, • ai][ ?- q_a,-a, - (c[-a,)' L ]- . 

. . - . J . . . _ j 2 . «, j xa, - aus, j fxal)^- 2-<2i f xai-5/jfS'j. )2x4,| 

| a, x « i i ( 

p ^ 

r — ^ ‘ ^ 2 *<2/ I ^ xa 2 ./ *t ^ 2* a 2 }^* &>| 


(2.21) 
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or 


o?j>_ _ itctr-n. ^ l- 

*>» dl |4 x a, l 3 - Tl K*^| L 

( 2 , 22 ) 

Finally, combining Eqs. (2.13) and (2.18) yields 


p • p 0 Id =. 

'D l? 


(l )'a, -> P^z + P 3 ^)- Volz - 

p 1 ® ^ jp- 1 2 |l .+ p ^ 2 xs ^ 

-a tf, 7 a t aaj 

T f - -Oft - - 7 Zt> — 
s) ■ j ~ 0.2.x n —=■ -t n xfl, - — — + at x a i 


OtZ 


Da; 


51a I 

7 7 ) J 


Or 


7 n J } a, xa^ 
(2.23) 


*3 7a _ P 


^ p |«t nail <■ 

"n < a , 




151 |ixCU 


/ 2’fli *«i 

. 1 * &{ \ 

' 151 

ll*5, (V 


+- 


(a'ta)[ «.■<*, . 4xdi-« 4-a t ] 

ia,xaijL —7-7 ^7=nr ■+ -7T] nr.Tr.iW 


5 i s- ji- 


lt I l4xfl,l 
7 -f -5, 


11 ^ ^ TT ^ *- 2. x / 

lixa 


Uf 

) 


I 4 X ^ J 1 i 111 


(2.24) 
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$ince, according to Eqs. (2.13) to (2.16), and Eq. (3.44) of 
Ref . 3 , 

a, li-dil 


K<*t\(- a . 1 *-a x +a. i q_-a > )z 


- - ^ (q'a a 3 ) 


I 


a' xa 2 . 5 1 , 


- q_xa z 


(2.25) 


and 


— := — (-n * a. f a- Qi kci z - a/XUzg-ajXrt)- 

\&t Ka^\ 

jOj xa z \ {- a z ci‘ a 2 -f a 1 £ . a 1 )- 


= -9_x (a’ L xa i ) -= - L _ = -£ 


^ ^ j 


a 1 ■ <a . * a 

(2.26) 

Equation (2.24) is equivalent to the desired expression 
for V 0 1 $ i 




— - <i_ < d 


1 


|4.x 


19. x ^,r / n( 


(2.27) 


■y , A ^Alternative Proof 

In order to verify Eq. (2.27), note that according to 
Eqs. (2.9) and (3.50) of Ref. 3 
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2* * * T- 

^ Zv 


= -2.7 


'ijvf 


-£[ r. (-0f) -- K ( 


^ 7^(t) 1<V<M 


(2.28) 


Noting that 
OJ 


/Zx a , ) ~ -p— ( P o ■* *} K J* ( P/ + 7 ^ 3 ) - 0 


~ (****)- ;£- ^o + lP;)/(P^fP3)^o 


7 


o 


f 


(2.29) 


yields 


2 ■ a* 2 * # 2 


/Z) 

■o 

Vf 


1 |?| | 4 *«il / 

f r *iJL + z.a t (-.JjLL-)l 

c[\ *• is I J / J 


lil 3 /J li^av ) 2 
1 (4.1 L J | q /a. I r 


^ 2 x 0.2 _ 


* . _ . , ~ -~~r (<2|X«i-f?Xp 3 )+ i * ^ 2 3 , 1 ’ . iV. ) ~ 

®f M 3 | 4! 2 1 )^|? j 

-[£•"£ f a,x?a * r*pjV 3 1x5:,. f* l { ]-2 


\H- 


(2.30) 


and, similarly, interchanging indices. 


9 


>7 V 


9 -- &i 


v |£/ l4*«J 


2 x ft; \ __ 

0 v n . I 2 * / 
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(V 1 i x a, + £* Pi) - % 


a, 



l 

\l \ 5 


(2.31) 


and thus 

ZfSlj-A. 2 -*'* 5 ~'*^ +3 rfxa.f.^- 2x^£-a,j]^p: 

(2.32) 


On the other hand 


i'- & 


= - U • F, 


< 2 . • « , x _ 


m 3 


u r ■£•«'*** 

F-[£<Z a, xa t -3 £-d, *di £ ] — — . - 
- - V-[_2- 4*? a,xaj.-t3 (5*a, £-5i - £ x a L £- a , )]— 


HI 5 


(2.33) 


5ince 

£•£ - £-7 £• a, * Tj. = 

2 X (''7, x 5V J - £* v 7 = 

£ • «2. 7, -fx *>- 1 •«,«*• £x7 - 
=- £ • a» 9 • \ x a , + £- a,y l - 9 .xa 2 . 


(2.34) 
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Since 9 is an arbitrary vector, Eq. (2.25) is equivalent to 


-v 0 ~ (- 1 J -ix*i 

(2.35) 

Equations (2.32) and (2.35) are the desired proof of the 
validity of Eq. (2.28). 


2.5 -Singularity at g = 0 

If the point P Q belongs to the surface 21^ , the integral 
in Eq. (1*7 ) is singular. In the following, the type of 
singularity is analyzed and it is shown that the principal 
value of the integral must be used. Consider a small circle 
of radius £* in the neighborhood of the singularity. Assume 
that the point Pq is at Very small distance from the surface 
and consider a small circular element 21 ^ on with 

the center on the normal projection of Pq on the surface z 4 
and radius £ . 

Assuming the z-axis to be directed along the normal n, 

Eq. ( ) reduces to 


K - *■ jj 5 * ^) d £ + «£ ■ 

with (for symmetry reasons, the derivatives 
Xq and yQ are zero) 



n i 


o 


(2.36) 

with respect to 


/ttv* TF '"" 


P 0 is not on the boundary of . 


O 
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* 2 '"’ *7 (t'f'iof 17 - ds 3 *7 ~ 


2// /i - 


' + *«■ I__ 

//v $7 ° 


- 


_ £ 


- ir^PTF 1 . 




<— •/— ^ 

--Zlin c 


(£ v ? 0 7 /z 


As Z n approaches zero, one obtains 


(2.37) 


If &(+)<** = Petd)** 


with 


(2.38) 


- p^(r¥r.-- ! 

^k _ (2.39) 

It may be noted that the first expression in Eq. (2.38) 
is not singular. Hence, Eqs . (2.8) and (2.27) (obtained by 

using the first integral in Eg. (2.38)3 are still valid even 
if the point is on the surface. 


2.6 General Element 

In this subsection, it is shown how the results obtained 
thus far can be rewritten in a more expressive fashion. For 
the sake of simplicity, introduce the following notations 
(Fig. 3> 

* P 0 is not on the boundary of YL& 
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£(0 ll ~ Ppp ~ P - <5 i 

l(-u) ~ P^p- p a) - Oi 

a (-».-»)- P r ^ = q 3 
*L (l 1 '0 - PpM ~ P C?4- 

(2.40) 

Note that 

a, fy-n) .(54- fli)/2 
<2 V " I } ~(<3 1 — £ 

Oa)/2 

(2.41) 

Next, combining Eqs. '(2.7), (2.27) (2.40), and (2.41), one 

obtains 

</*< = u<iJ)-iv (-h0+ T v (-o -i)-Tv (h-0~ 


<2. * # 2 . 

2 X 4% 

1 • «' i 

lx a 1 1 


Ill 

1 X i. 1 ^ 

/ 7T ! 
i i.( 

Cl ( | J 

I*'» T'i 



^ * # t 

I X a, 1 


it. i 

l^xaVl 1 

)4 \ 

1 ^ J 


r .z-**. 

I *5*. 

1 - a, 

7x5, 


L lii 


17i 

l ^xa,! 1 

J f '-hV 

f l-ai 

°i xfl I 

7-a, 

ixa, 

\ 

l \l\ 

l^xOCjt 

li f 

17*5:, t . 

r- 


(2.42) 
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or ^ 

Hf, t - -^Qi -&■»} gfXfg~,-gV) Qr(ti,-Gt) Q,x(oj-Qi) 

■ I3il | Q-oi |q,| la, xra,-«l)| l 

— &*■ ' 6?i: X. ^ gj^-<Sl. i) + Ql'(Q, -Qi) Qx K (&)-&>. ) 

l&xl \Q-l\ |5^x(ffi-5jl l 

_. a j • ( 5 3 . ~ <3 3 ) Qj * ((2x.-&l) fl 3'^<?4 -$ 3 ) 0; X ~Ql) 

IQj| i 5 3 1 

— ^4 C < 3 1 - $4 J <54 x ( Qi - Q4) . -Sj) GIa-x (Q4 -Q3) 

l a 4 | I^r fcUQ+'Qjf 


(2.43) 


or 


^hi = 


&±kQ 1 r 
0.4 * 0) i 2 
ffj X & 1 


Qa- * Qa - Q \' , &\ • (2 1 - Q\ • ©4 

■ — 1 T 


KS 4 .I 


Ifiil 


\ &■[ x ( 2 i I 2, 

51 x & j 

| G v * 

5*3 x ©4 


Q \ ■ Q i - 6! * ' Q i , r S a. ~ O r Qx. 
— — — - — — 


is. i_ 

r Qr^ • 6^ ^ Q 2 . ■ (A 3 


©3 X G*l z L 


\ © 2 . 

1Sl3\ 


+ 


J®i_ 

£•< < (S 3 2_° x ' a ^- 

\ ©2 i 


4 


-t 


+ 


$. 1 # ^5-3 ~~ ©3 * ©4 ^ Qla • Q4. — & 1 * <3 4 J 


| £M 

(2.44) 

It may be noted that each of the four terms depends upon two 
corners of one edge of the element. Hence, Eq. (2.44) is 


independent of the numbering used (it dependSj however, upon 
the direction of the numbering which is anticlockwise with 
respect to the normal n) . 
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Next, consider the limit of Eq. (2.44) when one edge 
shrinks to zero, that is when the hyperboloidal element 
reduces to a triangular element. As mentioned, the numbering 
is inessential. Hence, without loss of generality, it is 
assumed that 3 3 , (see Fig. 4-) . By setting 


(3 4 - $3 — A (3 4 3 

(2.45) 

where (3 4 ^ is a unit vector and h tends to zero. The last 
term of Eq.' (2.44) yields 

Q l k (^4 ( G 3 • ~ <3 3 ’ 34 ‘ - Q $ * & 4 \ _ 

*-■>0 i ¥ 3 * 1 an jbm. 

£iAM &3 * @4 3 _ /_ A <2 3 • 4 J 0 . 4 .- Ql 2 \ __ 

0 ?\ 1^3^ IQ 4 I 

6^-' > 0-3 1 *0.4.3) ' lC 2 i( i Q4A ' 

(2.46) 

Hence, for triangular elements 


*~h k 


cS i\ * Q\ f 

&3 ■ 

&3 - Qi ‘ & 3 t 

Qrflt-Q, • 


1 Q] t <3\'t* 


1'Qj 1 • 

I5,| 

- 

6^ x 0 2 ' 

■ ja_ 

\ * Ql 1 - (3 1 * & x 

c3 ~ <2i ■ 

•Or 1 

1 s-4* - 


^ . .. ._ .... .-. . — j_ 

I<3,1 


- 

(2iX 0-2 

i'A 


6 ? 3 • <3? — 0 . 


15** Qjf 

l 

1SJ 

1 •. 

„ 


(2.47) 
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Similarly, for a polygonal element with n corners 


+ 4, 3 -+ ' " ■+■ j n , 


(2.48) 


with 


7. - / g t --Qi -Qt'-fly QyQj-Qi'&j \ 

l/ lai-xSj ' 1 v |5- 1 ) 




I4YI 


(2.49) 

Equation (.2.48) can be proved as follows. The solid angle 
is an additive quantity. Hence, ^ , which is the gradient 

of the solid angle is an additive quantity. Thus, the 
general proof is obtained by mathematical induction: assumed 

to be true for n = n^, it is shown to be true for n = n^ ■ + 1, 
Thus (.see Fig. 5 for the case n^ = 4), noting that T t ‘j = - Ty 

T,,,. + +T 3,4^ 4,5 " ^4,1 


(2.50) 


in agreement with Eq. (2.48). 
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SECTION III 
SUMMARY AND RESULTS 


Consider Eq. (1.3); dropping the subscript 0 yields 

z B 

Assume that the geometry of the wake is prescribed as 
straight vortex lines or from the preceeding iterations 

(See Fig. 5; for a description of the iteration procedure 
see Ref . 5). Divide the wake into L strips, , each bounded 

by two streamlines. Divide the surface of the body into 
small polygonal elements (hyperboloidal quadrilateral, or 
triangular, for instance) . 

Then, Eq. (3.1) can be approximated by 


f- f »» II 

hZ ' Zk 


L 


0 


^ on 


(-)d'z. 


t 


(3.2) 


Next, assume that, in virtue of the Kutta condition it is 
possible to replace with the values of D at the centroid 

of the element having an edge in common with the strip Z^* 
Then, Eq. (3.2) can be rewritten as 


where 




(3.3) 


(3.4) 
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if 51^ has no edge in contact in contact with the wake, while 

+ (3.5) 

if has an edge in contact with the strip U- 

The perturbation velocity, 'iTjv > the centroid, fyi 
of the element is given by 


N 




where 


h / 


***=(>// £ ^ 


p= Pi, 


(3.6) 


(3.7) 


0 Jj dn <■ r 

Finally, imposing the boundary condition at the centroid 


of the elements, yields the system 


where 


while 


['M] I M " I B *' 


f\ h $ “ ^ k 


r^jL 

Un JP-Ph 


(3.8) 


(3.9) 


(3.10) 


is prescribed from the boundary conditions. Solving Eq. (3.8) 
yields the coefficient D^?: then, it is possible to evaluate 

/^through Eq. (3.6). 

The integral in Eq. (3.4) can be evaluated by using 
Eq. (2.48) for a general polygonal element^ or Eq. (2.4*3 for 
triangular elements, or Eq. (2.44) for hyperboloidal quadri- 

Note that, if the element includes 


lateral elements . 
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a strip it will be convenient to approximate it with 

a series of quadrilateral subelements. Then 21/^ can be 
treated as a single polygonal element: in this way the con- 

tribution of the edges (which would eventually eliminate 
each other) need not be evaluated. 

This formulation has been implemented into a computer 
program, ILSA , (acronym for Incompressible Lifting Surface 
Aerodynamics). See also Ref. 5. Figure 7 shows the lift 
coefficient distribution per unit angle of attack for a 
rectangular wing of AR = 8 , at Mach Number M = 0. A conver- 
gence study for various numbers of wing elements is also 
shown and compared to the result obtained by Yates (Ref. 6) . 

The results obtained with ILSA indicate good agreement with 
existing ones and a fast rate of convergence. As mentioned 
before, a better wake geometry can be obtained by an iteration process. 
This process is shown in detail in Reference 5. 
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